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Abstract 

In this chapter some results related to Shunt Active Filters (SAFs) and obtained by the authors and some 
coauthors are reported. SAFs are complex power electronics equipments adopted to compensate for cur- 
rent harmonic pollution in electric mains, due to nonlinear loads. By using a proper "floating" capacitor 
as energy reservoir, the SAF purpose is to inject in the line grid currents canceling the polluting har- 
monics. Control algorithms play a key role for such devices and, in general, in many power electronics 
applications. Moreover, systems theory is crucial, since it is the mathematical tool that enables a deep 
understanding of the involved dynamics of such systems, allowing a correct dimensioning, beside an 
effective control. As a matter of facts, current injection objective can be straightforwardly formulated 
as an output tracking control problem. In this fashion, the structural and insidious marginally-stable 
internal/zero dynamics of SAFs can be immediately highlighted and characterized in terms of sizing 
and control issues. 

For what concerns the control design strictly, time-scale separation among output and internal dynamics 
can be effectively exploited to split the control design in different stages that can be later aggregated, 
by using singular perturbation analysis. In addition, for robust asymptotic output tracking the Internal 
Model Principle is adopted. 

In authors' opinion, SAF case is an illustrative example of common issues in dimensioning and control 
of complex power electronics equipments, hence the proposed design approach can be generalized for 
such class of systems (e.g. large-power electric drives; generators and converters, particularly for re- 
newable energies; Uninterruptible Power Supplies, UPS, and power supplies for special applications as 
particle accelerators). Remarkable the role of "system theory approach" in enlightening crucial sizing 
issues, this fact strongly testifies how relevant the "control viewpoint" is in all the fields of engineering, 
particularly when complex dynamic behavior is requested. 

The presented chapter has been invited for possible publication (after a review process) in the book 
"Robust Control / Book 1", ISBN 978-953-307-421-4, by INTECH (www.intechweb.org). 

1. Introduction 

Harmonic pollution in the AC mains determines additional power losses and may cause malfunction- 
ing or even damage to connected equipments. Distortion of the currents circulating on electric mains 
is mainly originated by non linear loads, as AC/DC uncontrolled rectifiers used for motor drives, that 
absorb undesired current harmonics. Therefore, local countermeasures have to be taken in order to keep 



the portion of grid affected by distortion as small as possible, hence preventing relevant power losses 
and "saving" other equipments, connected to the rest of the grid. 

Traditionally, passive filtering components have been adopted to cope with harmonic compensations, 
however they are affected by several drawbacks; they are very sensitive to network impedance varia- 
tion and environmental conditions, moreover they need to be tuned on fixed frequencies. In order to 
overcome those limitations, in the last decades, thanks also to the fast growth in power electronics and 
control processor technologies, a remarkable research attempt has been devoted to the study of the so- 
called Active Power Filters (APFs), both from a theoretical and technological point of view (see Gyugyi 
|& Stryculal ( |1976) , [Akagil ( |1996| >, |Singh & Al-Haddad|fl999) ). These devices are able to properly work 
in a wide range of operating conditions, providing better performance and overtaking intrinsic limita- 
tions of passive devices, they are far more insensitive to network impedance, they can be tuned onto 
different frequencies just varying some software parameters. Furthermore, the system reliability is im- 
proved, resonance phenomena are avoided and a diagnosis system can be implemented on the control 
processor to monitor the system variables and adopt some recovery strategy in case of faulty conditions. 
In this chapter, the general issues related to analysis, dimensioning and control of a particular class of 
APFs, the so-called Shunt Active Filters (SAFs), are addressed; the main purpose of this kind of power 
system is to inject into mains a proper current, in order to cancel out, partially or totally, the power 
distortions generated by nonlinear loads. The SAFs considered in this work are based on a three-phases 
three- wires AC/DC boost converter topology (see Fig. [TJ connected in parallel to the distorting loads. 
The first step to properly design a SAF is the selection of suitable hardware components; as it will 
become clear in the next section, owing to the structural properties of the system, the sizing procedure 
cannot be considered apart from the canonical control aspects, hence a correct dimensioning algorithm 
jRonchi &~ Tilli 2002) is proposed to ensure feasibility of the desired control objectives and to minimize 
costs. In addition, according to such method, it is shown how a time-scale separation between different 
dynamics of SAF usually takes place "for free". This point is very useful for control design and stability 
analysis. 

Once a correct hardware sizing has been carried out, the first control issues to deal with are: the cur- 
rent/power control algorithm and the load current analysis method adopted to define the filter current 
reference. Various solutions have been proposed in literature. As regards current/power harmonic 
tracking, in ^Chandra et al. 2000[> an hysteresis current control {Kazmierkowski & Malesam] |1998| > 
is proposed, while in ^Jeong & Woo| |1997} predictive current control is adopted. For what concerns 
the generation of the filter currents reference, beside Fast Fourier Transform Techniques, instantaneous 
power theory ( Akagi et al. 1984 1, time domain correlation techniques ( [Van Harmelen & Enselin|["l993! >, 
notch filter theory ( Rastogi et al. 1995| and other methods have been proposed. Solution based on state 
observer have been proposed, too, as in (Bhattacharya et al.|[T995> and frilli et al.|[2002| >. 
However, what renders the SAF control problem challenging and different from other conventional 
tracking problems is the presence of peculiar and unstable internal dynamics, given by the voltage dy- 
namics of the DC-link capacitor bank. This capacitor bank is the main energy storage element, which 
provides the voltage, modulated by the control, to steer the filter currents and, at the same time, is re- 
quired to oscillate to exchange energy with the line and the load to compensate for current harmonics. 
Actually, this element needs to be carefully considered also in the previously-mentioned dimensioning 
stage; a correct capacitor sizing is crucial for control objective feasibility, whatever control technique 
is adopted. Moreover, it can be shown that, if perfect harmonic compensation is achieved, the DC-link 
voltage dynamics are unstable, due to the system parasitic resistances that lead to a slow discharge of 
the capacitor. Hence, a suitable stabilizing action for DC-link voltage dynamics needs to be provided. 
Since no additional circuit is used to feed the DC-link capacitor independently of the three-phase port 
used to inject currents, (see Fig. [TJ, the voltage stabilization would need to be integrated with the con- 
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Fig. 1. Shunt Active Filter scheme. 

trailer devoted to harmonic compensation (the AC/DC boost-based SAF is an underactuated system). 
This is a crucial point and it has to be tackled preserving harmonic compensation performances as far 
as possible. 

In this work a power/current controller, based on Internal Model Principle, (see Mar coni et al.| (2003 ), 
Marc oni et al.| (2004 1, Mar coni et al.| |2007 1) is designed in order to cancel current harmonics, ensuring 
robustness with respect to SAF parameter uncertainties. By exploiting the internal model approach, the 
proposed solution also allows to merge and solve at the same time the two above-mentioned problems 
of current harmonics isolation and current reference tracking. As regards the robust stabilization of 
the DC-link voltage internal dynamics, a cascade control structure is proposed. An additional volt- 
age controller, acting on the references of the power/current controller, is introduced. This controller 
is designed taking into account the structural voltage oscillations required for harmonic compensation 
and minimizing the impact on harmonic compensation. In particular, by exploiting a proper averaging 
(Sanders et al. 1991) of the capacitor voltage dynamics, the average value of the capacitor voltage is 
chosen as output variable to be controlled (Hanschke et al. 2 006| >. 

As far as the overall stability is concerned, the previously mentioned time-scale separation between 
portions of SAF dynamics can be effectively exploited to decouple power/current tracking and voltage 
stabilization control problems, using averaging and singular perturbation theory techniques (Khalil 
[1996) . 

This chapter is organized as follows. In Section|2] the general framework is described, the SAF model 
is derived and the control objectives are formally defined. In Section[3] two methodological approaches 
are presented for the SAF components sizing. The first one is based on the knowledge of the load cur- 
rents harmonic spectrum, the values selected for the hardware components are the minimums allowing 
the SAF to deal with the considered load distortion. Differently, the second approach is related to the 
maximum current of the AC/DC boost switching devices. In this case the selected components values 



are the minimums which enable the SAF to compensate for all possible loads giving distorted currents 
smaller or equal to the switches peak value. In Section[4]both the internal model-based power/current 
controller and the averaging voltage controller design are presented, stability analysis is carried out re- 
lying upon the time-scale separation imposed by the design algorithm; both the power and the voltage 
subsystem are proven to be asymptotically stable, then practical stability of the overall system is claimed 
exploiting general results on two time-scale averaged systems (Teel et al. , 2003). The effectiveness of 
the proposed control solution is tested in Section|5]through simulations. 

2. Shunt active filter model and control problem statement 

The scheme of the shunt active filter considered in this chapter is reported in Fig. [T| as mentioned in 
the introduction it is based on a three-phase three-wire AC/DC boost converter, where the main energy 
storage element is a DC-bus capacitor, while the inductances are exploited to steer the filter currents by 
means of the converter voltages. The switching devices of the three-leg bridge (also called '"inverter"') 
are usually realized by IGBTs (Insulated Gate Bipolar Transistors) and free-wheeling diodes. 
In this work the following notation is used to denote the SAF variables; v ma b c ={y ma , v m b, v mc ) T is the 
mains voltage sinusoidal balanced and equilibrated tern, i m =(imaJmbJmc) T are the mains currents, 
h = {ha'hb>hc) are me l° a d currents, while i=(i a ,ii,,i c ) T are the filter currents. L indicates the value of 
the inductances, and C the DC-link bus capacitor value. 



2.1 Mathematical model 

Considering the inductors dynamics, the filter model can be expressed as 
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where R is the parasitic resistance related to the inductance L and to the cables, vnk is the voltage 
between the nodes N and K reported in Fig. [T] v(f) is the voltage on the DC-link capacitor, and 
u \ = {u x ,Uy,u z ) T is the switch command vector for the legs of the converter. Since a PWM (Pulse 
Width Modulation) strategy is assumed to control the inverter, the above-mentioned control inputs can 
be considered such that «i; £ [0, 1], i = X,y,z. According to the three-wire topology for any generic 
voltage/current vector x it holds 

£ x, = (2) 

i=a,b,c 

hence, from the sum of the scalar equations in |TJ it follows that 

u x (t) + u y (t)+u z (t) 
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it can be verified by direct computations that 
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[1 1 l]u abc (t)=0 Vf>0. 



(5) 



For what concerns the state equation relative to the capacitor voltage dynamics, it can be derived con- 
sidering an ideal inverter and applying a power balance condition between the input and the output of 
the filter, then replacing {3} into {T}, the complete filter model results 
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where the filter currents dynamics have been written in a more compact form with respect to QJ, mul- 
tiplying the current vector by the identity matrix of suitable dimension /}. Exploiting equations 13), 
j5), the system model can be reduced to the standard two-phase planar representation of a three-phase 
balanced systems ( Krause et al. 1995 1, which can be obtained applying the following coordinates trans- 
formation 
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the SAF dynamics expressed in this a — /3 reference frame become 
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according to the hypothesis of three-phase balanced sinusoidal line, the ideal main voltage tern can be 
expressed in the above-defined bi-dimensional reference frame as follows 

[v ma v„,p] T = V m [cos{(O m ) sin((O m )] T 

where V m is the voltage amplitude and co m the grid angular frequency. For what concerns the control 
vector u abc , in this reference frame the eight possible configurations of the switching network (reported 
m Tab. [[} can be mapped in the a — j8 plane, obtaining the vertexes and the origin of the feasibility 
space illustrated in Fig|2] while each point in the hexagon can be obtained as mean value in a PWM 
period. As it will become clear in the next sections, in order to simplify the control objectives definition 
and the controller design, it is very useful to adopt a further transformation from the two-phase current 
variables [i a ip ] T to a two-phase real-virtual (imaginary) power variables defined as 
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where 
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In this so-called synchronous coordinate setting, aligned with the mains voltage vector, the model of the 
SAF is expressed as 
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Table 1. Control function table. 




where 

3 2d 
. e = 7^ — . ^mrf = V m , = q T a pu a p (11) 

it is further to notice that, since the filter currents, the mains voltage and the DC-link voltage are mea- 
surable, the full state (x,v) is available for feedback, moreover the actual control action u = [u x u y u z ] 
can be determined from u a i, c , which in turn can be derived from 

As regards the load description, the same two-phase real- virtual power representation can be used, in 
particular following (Aka gi et al.||1984| >, the load currents can be approximated as periodic signals 
given by the sum of a finite number N of harmonics, with frequencies multiple of f m = a>,„ / 2%. Hence 
the load currents can be expressed in power variables as 

N+\ 

x lj = X lfl) + Xljn cos ( nC Omt + Vjn), j = d,q (12) 

«=1 

where the harmonics amplitudes X^q, Xj„o, 3Qdn> Xi q „ and phases y/^„, y qn are constants. Since the load 
currents and the mains voltages are measurable, also the variables (x/j, x\„) will be considered known 
and available for control purpose. 
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2.2 Problem statement and control objectives 

Roughly speaking the main control objective of the considered SAF is to steer the variables xj, x q , 
injecting power into the line to compensate for the load harmonics. However the ability of tracking 
current references relies upon the energy stored in the DC-link capacitor, which is the main power 
source of the filter, therefore another general objective is to keep the DC-link voltage confined in a 
suitable region, to avoid overcharge and, at the same time, to ensure the capability to steer the filter 
currents. On the other hand the ability of maintaining DC-link voltage into a suitable region is strictly 
related to the power exchanged with the mains, which in turn is affected by the current harmonics to 
be compensated for. The general control objective is then two-folds; one related to the tracking of 
current disturbances, the other concerns the voltage internal dynamics stabilization. In this paragraph 
a precise and feasible control problem is formally defined, recalling the considerations made above, 
and assuming that a suitable dimensioning, that will be deeply discussed in the next section, has been 
carried out. 

Bearing in mind the power variables representation of a generic nonlinear load expressed in l |12[ >, it 
turns out that the only desired load component is X^o, since it represents first-order harmonics aligned 
with the mains voltages, while the remaining part of the real component x^ — X^o is an oscillatory 
signal with null balance over a line period, and the imaginary component x\ q represent a measure of 
the misalignment between mains ideal voltage and load currents (see Mohan et al. ( 19891) and do not 
contribute to the power flow. In this respect, the terms x/j —Xmq, x\ q are undesired components which 
should be canceled by the injected filter currents, hence ideally the control problem can be formulated 
as a state tracking problem, for system < | 1 0| >, of the following reference 



a prefect tracking of this reference would ensure pure sinusoidal mains currents perfectly aligned with 
the mains voltages. However, this ideal objective is in contrast with the requirement to have a DC-link 
voltage bounded behavior. In order to formally motivate this claim, consider the steady state voltage 
dynamics in case perfect tracking of the power reference x* (t) is achieved, after some computations it 
results 



the signal eY(x* (r)) which drives the integrator is periodic with period T = 1 //,„, and it is composed 
by the sum of a zero mean value signal eL^cIq —x*) t x*, and the signal eL(M(R,L)x*) t x* which has 
negative mean value as long as parasitic resistance R or reference x* are not zero. By this, no matter 
the starting voltage value of the DC-link, the capacitor will be discharged and the voltage will drop, 
providing a loss of controllability of the system. 

To avoid this phenomenon, the reference must be revised, taking into account an additional power 
term, which should be drained from the line grid by the active filter, in order to compensate for its 
power losses. Following this motivation, and recalling that the unique useful component for the energy 
exchange is the real part of the power variables, the ideal reference signal l |13| l is modified as 
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in which <po is a solution of the following equation 




(16) 



this represents the power balancing condition which guarantees that the internal voltage dynamics in 
case of perfect tracking of the modified reference x^ is 

^=eT(4,(0) (IV) 

and the right hand side £ x ¥(x^ )g ) is periodic with period 1 / f m with zero mean value. A brief discussion 
is needed for the solutions of equation \16\ , it has two real positive solutions if the following condition 
is verified 

E} nd > *R Z fmf o M +X* q 2 (z))dZ (18) 

from a physical viewpoint relation fTSt set an upper bound on the admissible undesired components 
which can be compensated and on the parasitic resistance R, however, as typically E m( j » R, this 
condition is not limitative at all. The two solutions of {16} under condition l | 1 8[ i are 
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the first solution, minimizing the power drained from the line grid to compensate the power losses, is the 
physically most plausible, because the power consumed by parasitic resistances in the filter is usually 
quite small, hence it will be considered throughout the chapter. 

The control problem which will be faced in this work can now be precisely stated; the issue is to design 
the control vector u a i, c in a way such that the following objectives are fulfilled: 



A) Given the reference signal xJL defined in 1 15 i, asymptotic tracking must be achieved, that is 



lim(x(f)-x^)=0; (20) 

B) Given a safe voltage range [v m , vjvr], with vm > v m > 0, and assuming v(fn) £ [v m , vm], it is required 
that 

v(f) € [v m , v u ],Vt>to; (21) 

it can be verified that the tracking of the modified power reference is potentially achievable keeping the 
voltage dynamics inside the safe region, only if the zero mean value oscillating component of eY(jc^) 
is properly bounded, this can be ensured by a suitable capacitor design. 

In the regulator design, saturation of the actual input u\, imposed by PWM strategy, will not be taken 
explicitly into account, also this approximation takes advantage of a correct sizing methodology; as it 
will become clear in Section|3] a suitable choice of the DC-link voltage lower bound v m , depending on 
the currents to be compensated for, has to be made to meet the constraint u\ £ [0, 1], at least when the 
power tracking error is reasonably small. 

A further consideration needs to be made on the requirement v(fo) £ [v m , vm] ; according to the AC/DC 
boost converter theory (Mohan et al. , 1989 1, the natural response of the system would lead the DC-link 
voltage at twice the line voltage peak value, due to the resonant behavior of the LC pair and the free- 
wheeling diodes of the switching bridge. If a proper design has been performed, this value is expected 
to be greater than the voltage range lower bound v,„; hence, after a transient period, the controller can 
be switched on having the initial voltage value inside the admissible region as required by objective B. 



Finally it is further to remark that x^ depends on parasitic resistance R through (16, hence it has to be 
considered as an unknown variable, to be reconstructed by estimating the power losses by means of a 
proper elaboration of the DC-link voltage signal. 



3. Shunt active filter sizing methodology 



The aim of this section is to present a precise algorithm to properly select the SAF hardware compo- 
nents, two different design objectives are considered, the first is to select the minimal component values 
dependent on the level of current distortion imposed by the load, while the second is to find the min- 
imum capacitor value necessary to compensate all the possible loads compatible with the maximum 
current rating of the inverter switching devices. Both the methods are control-oriented, that is they en- 
sure the feasibility of control objectives stated in |2.2| and that control input saturation is avoided under 
nominal load and line voltage conditions. 

The proposed design method is based on the model derived in Section|2] a further approximation is con- 
sidered with respect to equation d6); the inductors are modeled as pure inductance, that is the parasitic 
resistance R is neglected, while ideal mains voltage tern and converter switches are considered as in the 
previous section. 

3.1 Inductance value selection 

The inductance value can be selected regardless the loads, hence this part of the design procedure is the 
same for both the design objectives previously defined. 

The design criterion is based on the maximum current ripple AIm pp allowed for the filter currents; 
current ripple is a consequence of the PWM technique applied to obtain the reference command value 
u *abc' ^ nas to b e bounded in order to limit high frequency distortion. The actual command vector 
u abc{t) an d filter current i(t) are affected by a ripple component 

((f) = f(t) +Ai(t) 

(22) 

Uabc(t) = u* ahc (t) + Au abc (t) 



substituting these expressions in the state equation d6l it turns out 



L — — = -Au„ 



dt -LlabA'Ht) (23) 



by simple computation it can be showed that the worst ripple case occurs when the desired command 
value u* bc is in the middle of a feasibility hexagon side (see Fig. In this condition, assuming that 
the DC-link voltage has constant value V in a PWM period, the peak to peak current ripple is 

AI f' +Ts/2 dAi(t) V 24 

PP Jt dt 6/pwmL 



where the sampling period T s and the PWM frequency fpwM are assumed already set before starting 
the sizing procedure. If the peak to peak ripple must be bounded by the desired maximum value AIm pp , 
the following inequality needs to be fulfilled 

L > ^ => L min = V »- (25) 

6fpWMAlMpp 6fpwM^Mpp 

the upper bound of the voltage range vu depends only on the kind of capacitor and it can be sup- 
posed already chosen before starting the design procedure, hence the minimum inductance value Lmin 
compatible with the desired maximum current ripple can be selected applying equation <[25). 
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3.2 Load-based approach 

Let us now consider the first design algorithm based on the knowledge of the load to be compensated 
for. The load distortion will be modeled as in equation (12) , taking into account the constraint on the 
maximum current I max of the device implementing the bridge switches. The switching devices sizing 
depends on the total amount of power (distorted and reactive) P = 3V„,rmsIsafrms mat me filter has to 
compensate for (if the load is known then P is known), hence by the route mean square value Isafrms> 
the maximum current that the switches need to drain can be readily obtained as I max = V2Isafrms- 
The desired filter currents (denoted with *) necessary to fulfill the tracking objective A defined in |2.2| 
can be effectively imposed by the converter if each component is less than the maximum allowed value, 
i.e 

?(t) = [«*(0 e h {t) f c {t)] T < w[i i i] T , Vf (26) 

this feasibility condition can be graphically represented considering that each projection of the filter 
currents vector must be less then l max , hence the feasibility space is an hexagon similar to that reported 
in Fig. |4](obtained taking P = 45kVAR as filter size, V,„rms = 220V and then I max = 10A). Therefore 
condition {26} can be readily checked considering the inscribed circle in the feasibility hexagon. If the 
load currents do not satisfy constraint l |26[ > the number of current harmonics to be compensated for has 
to be reduced, differently, when the filter performance cannot be decreased, the opportunity to connect 
two shunt active filters to the same load can be considered. 

Assuming that an inductance value such that L > L m ,„ has been selected, the voltages at the input of the 





six switches bridge can be calculated as 
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the above equation is obtained by inversion of equation {TO} with R = and expressing the model in 
the synchronous reference frame in current rather than in power variables, in order to directly consider 
the load currents in the design approach. 

The constraints on the command inputs need to be considered too, by l |27[ > the inductance value must 
be as low as possible in order to make u* dq feasible, taking into account also the current ripple limitation 
we select L = L m \ n . As mentioned, the choice of the of the capacitor voltage lower bound value plays 
a key role to avoid saturation issues on command inputs, this can be easily verified approximating the 
hexagon in Fig. [2]with the inscribed circle. In order to avoid control action saturation (assuming perfect 
power tracking) it must be imposed that 

,,„* I, _ \KbM\ „ Kbci*) 
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with ri n the radius of the inscribed circle and v* bc = v(t)u* bc . From 1 28 l design equation for v m can be 
obtained 
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usually v m is oversized with respect the value given by the inequality above, in order to avoid saturation 
even if non-zero tracking errors are present. If condition {29} cannot be satisfied, some alternatives 
need to be considered; the capacitor can be changed in order to adopt an higher upper bound vm, when 
the costs of the project have to be limited and the kind of capacitor cannot be substituted, the number 
of harmonics considered must be reduced until {29} is satisfied. To preserve the number of harmonics 
to compensate, the inductance value can be reduced, penalizing the current ripple and then tolerating a 
greater high frequency distortion. 

The capacitor value can then be selected assuming an ideal converter and writing the balance equation 
between the instantaneous reference power at the input of the six switches bridge and the power of the 
DC-link capacitor, hence 



(30) 



Pfilt(t) = [Vdq{t)] idq{< 

the corresponding energy can be calculated as 

Efut(t) = [ Pfm{t) 

by the hypothesis of sinusoidal load currents and ideal mains voltages Em t (t) is periodic of frequency 
f m and its mean value is zero. Defining 
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Vref=^— 

and imposing that the voltage variation corresponding to E max is v re f 
equation can be written as 

2-£mflA' 



c = 



ref 



(32) 

the capacitor value design 
(33) 



3.3 Switches-based approach 

As stated at the beginning of this section, a different design method aims to find the capacitor value that 
makes the filter able to compensate for the worst load compatible with the switches maximum current. 
If the resulting capacitor value is not too expensive, this method allows to design the filter only knowing 
the amount of current that has to be compensated. 

During the optimization procedure the load currents need to be the only varying parameters while all the 
other values must be fixed. The inductance value is chosen equal to the minimum compatible with the 
allowed ripple, while the minimum capacitor voltage v m is supposed sufficiently low to make simple the 
voltage control, and, at the same time, the resulting capacitor value feasible. Writing the filter currents 
spectrum in the d — q synchronous reference frame, an expression similar to \12\ can be obtained 

N+l 

0(0 =I j0 + X I JnCos(2mf m t + rirj n ), j = d,q (34) 

the parameters to be varied in order to calculate the worst E max are the (22V + 1) + 1 magnitudes and 
the 2N + 1 phases, so the following optimization problem 



nax| / [v dq (z)] T i dq (i;)di;\ 



(35) 



has to be solved with respect to the array z of 4(JV + 1) +1 variables, taking into account the following 
constraints 

• switches currents must be less than the maximum allowed, that is the current vector must be 
inside an hexagon similar to that reported in Fig. [4] This can be easily checked approximating 
the hexagon with its inscribed circle; 

• the control output must be feasible, that is the vector u a i, c must be inside the hexagon reported 
in Fig. [2] This can be easily checked approximating the hexagon with its inscribed circle; 

• harmonics components phases have to be greater than — % and less than K. 



Once E™' x has been determined, substituting its value in 1 33 i, the capacitor value relative to the 



switches-based design approach can be selected. 

In the discussion above, ideal mains voltages have been assumed, if also the grid line voltages are 
distorted, the capacitor has to provide more energy to the load, hence its value must be higher than the 
one calculated under ideal conditions. 

In case of ideal mains voltages the load instant power is the one calculated in l |12} and the only power 
term that the filter must deliver is x;^„ = E^Zj Xi ( i„cos(nCQ m t + ydn) — Vmiidn having zero mean value. 
If the mains voltages are distorted, their representation in the synchronous reference frame is 

Vmdq(t) = [Vm + Vmdn* VmqY (36) 



line voltages harmonic perturbation produces additive terms in the load instantaneous power expression, 
that by direct computation can be written as 

Placid = v mdn{ t )ildn{t) + Vmqil q {t) (37) 



the above equation shows that the filter has to provide more power to the load, furthermore the power 
mean value in a PWM period can be different from zero. Hence also assuming that the mean value 
becomes zero in a finite time, the capacitor must be oversized with respect to the ideal situation, in 
order to accumulate more energy. 
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Fig. 5. Controller structure. 



4. Robust controller design 

In this section the control problem defined in |2.2| is addressed, relying upon a suitable capacitor value 
given by the procedure described in the previous section, the two interlaced objectives A and B defined 
in |2.2| can be approached individually by exploiting the principle of singular perturbation. Two inde- 
pendent controllers (reported in the block diagram of Fig. |5j will be designed. An internal model-based 
controller (IMC) is proposed in order to deal with the problem of robust reference tracking (defined in 
objective A) for the fast subsystems composed by the power variables dynamics, while an independent 
voltage controller for the slow DC-link voltage subsystem is designed to produce a reference modifi- 
cation r/ which compensate the unknown power losses term <po, allowing to achieve objective B. The 
averaged voltage value is chosen as the controlled variable, and a phasor variables representation is 
exploited to design the regulator, this choice allows for the necessary voltage oscillation during nom- 
inal operation, and improves the voltage dynamics behavior with respect to other proposed solutions 
jMarconi et al.| [2007). Stability analysis is carried out in two steps; the reduced averaged dynamics, 
obtained replacing the steady state of the fast subsystem into the slow voltage dynamics and carrying 
out the average value to obtain a phasor variables representation, and the boundary layer system, ob- 
tained considering the SAF currents dynamics and an ideal energy storage element, are proved to be 
asymptotically stabilized by the proposed controllers. Then practical stability for the overall closed- 
loop error system is stated exploiting well-established singular perturbation and two time-scale systems 
theory results. 

Before detailing the proposed control structure, consider the first preliminary control law 

u(t)=v{t)u dq (t) (38) 



which is always well defined provided that v(t) > v m > for all t > according to objective B. Replac- 
ing l |38l > into |T0| > yields 

x = M(R,L)x+ -u + d 



dv _ T 

— — =£11 I 

dt 



now consider the modified power reference 



x* =x* + (T] 0) t (40) 

and define the change of variables 

x = x-x*, z = v 2 -V* 2 (41) 

where V* 2 = (v 2 t + vh) / (2) is the reference value for the square DC-link voltage. Note that the 
requirement B of having v(t) £ [v„, \>m] for all t > to can be equivalently formulated in the error 
variable z requiring z(t) 6 [— /* /*] for all t > to, with I* — (v|^ — v 2 n )/ 2. The complete system l |39[ > 
can be then expressed in the error variables defined in (41) , the transformed model results 

x = M(R,L)x — -u + do —x* +Mxt 

L n 1 (42) 

Z = £U T [x+xt]. 

The controller design will be carried out considering the error dynamics in l |42[ (, in summary the idea 
is to steer the closed loop dynamics toward a steady state in which z is free to oscillate within the 
admissible region, but its mean value is steered to zero (i.e the DC-link voltage mean value is steered to 
V* ), and x is steered to zero, i.e the power x follows a reference which is the sum of the term x* , which 
takes into account the undesired harmonic load components, and a constant bias Tj which is needed 
in order to compensate the power losses and to make the range [v m vj|f ] an invariant subspace for the 
voltage dynamics. 



4.1 Averaging voltage controller 

In order to fulfill objective B the voltage dynamics need to be stabilized, in this respect the subsystem 
composed by the capacitor voltage dynamics will be considered, a suitable reduced averaged system 
will be sought, and then a controller for the capacitor voltage DC component will be designed. 
The first step is to average the voltage differential equation to obtain the dynamics in the so-called 
phasor-variables, then, a control law, itself expressed on phasor representation, can be designed fol- 
lowing an approach similar to that proposed in ( Valderrama et. al 200 1 1, however in this work the only 
voltage subsystem is controlled using phasor variables, while the power subsystem is controlled in the 
real time domain. 

The controlled variable is chosen to be the time-window averaged value z a of the square voltage error 



Z, and the averaging is performed over the time interval [t — T,t\. In terms of i Sanders et al. 199 lp this 
average value is a zero-order phasor defined as 



z a (t) = [ z(r)dt 
Jr-T 



(43) 



the fact that z a is a zero-order phasor allows to obtain its derivative by simply applying the same aver- 
aging procedure to its differential equation in J42p 



1 ft 



f ~z(r)d 

Jt-T 



Za= h \ z(r)dT = el n T [x + xt ]dx (44) 



T Jt-T Jt-T 1 

note that the average voltage derivative can also be expressed as the difference over one period of the 
actual voltage, hence 

d . , d [' , s z(t)-z(t-T) 



this insight connotes the availability of z a for measurement in real time, and, as it will later clarified, it 
is of crucial importance for an actual implementation of the controller. 

All further elaborations will focus on the integral-differential equation |44j representing the averaged 
error voltage dynamics. This equation depends on u which is actually provided by the power tracking 
controller, to eliminate u consider that the differential equation for x in j42} can be rewritten as 

u = L(M(R,L)x* rl -x* +M(R,L)x + do-x) (46) 

replacing ^46\ into ^44\ the following equation is obtained 

Za = Y / T ( M ( R ' L ) x *r, -^+d ) T x^dt + eLD(x) (47) 

where £>(x) collects all the terms depending on the power tracking error x. The next step is to exploit 
the fact that the reference term x* is r -periodic (T = 1/ f m ), hence it results in a constant value when 
averaged over this period, this is a key advantage of the averaging approach for the voltage system. The 
T -periodic terms in j47} can be summarized to 

1 f' 

D* = - \(M(R,L)x* -x +d ) T x*]dz (48) 
T Jt-T 

since x* is periodic in T, D* is a constant disturbance, and, due to power losses induced by the parasitic 
resistance R, it also follows that D* < 0. For further simplification the integral operator can be applied 
to the occurring derivative terms. Using definitions {TTJ, l |40} , after some computations the averaged 
error voltage dynamics can be expressed completely in phasor variables 

L = £ [E md r\ a - 2Rv a - Lv a + LD* + LD] (49) 

where the following nonlinear term has been defined 

which enters l |49[ > with its average and its averaged derivative 

1 r> 

v a (t) = - jf v(r)dr 



(51) 



the averaged error voltage system is thus controlled by means of the averaged control input 

1.(0 = = f 1 ( 52 ) 

T Jt-T 

According to singular perturbation theory, the voltage controller design can be carried out considering 
only the reduced dynamics, obtained confusing the value of x with its steady state value x = 0. 
As previously remarked, this approximation can be introduced thanks to the small value of e which, 
multiplying the voltage dynamics in the second of \42\ , makes the voltage subsystem much slower with 
respect to the power dynamics in the first of \42) (this phenomenon is usually referred as two time-scale 
system behavior) that will approach the steady state much faster then z. Thus reduced voltage dynamics 



can be obtained by 1 49 1 simply dropping the coupling term D, because by definition D(0) = 0. 



The nonlinear terms v fl , and v a cannot be managed easily, beside non-linearity they contain an integral, 
a time delay and a time- varying term x* d . In order to simplify the mathematical treatment, a sort of 
linearized version of system l |49| > will be considered. This linear approximation is motivated by several 
facts; since the parasitic resistance R and the filter inductance value L are usually very small with 
respect to the term E md in every realistic setup, nonlinear term are much smaller than the linear ones. 
Furthermore the component x*, has no influence at all in averaging terms if T] is constant, thanks to the 
fact that it is T-periodic with zero mean value. Hence it will influence the averaged system only while 
tj is varying, and also in this case its oscillatory part will be filtered by the averaging procedure. As a 
result of the previous steps and considerations, the linearized averaged model for the reduced dynamics 
can be written as 

z a = £E md - <po] (53) 
where, as mentioned, q>o is the smallest solution of equation \16) . 

Now it is possible to design the control input r\ a in order to stabilize the origin of system {53), a standard 
PI regulator in the averaged variables is proposed 

r\a = -KpZa + 
. (54) 
9 = -EKrZa 

it is further to notice that the e factor in the integral action of the controller is introduced to keep the 
voltage controller speed in scale with the voltage subsystem to control, thus maintaining the two-time 
scale behavior of the overall system. 

In order to prove the asymptotic stability of the closed-loop system resulting by the interconnection of 
{54) and {53) consider the change of coordinates 9 = 9 — <po, which results in the closed-loop error 
dynamics 





Za 




9 



(55) 



£_ Za -E md K P E md 

dt [e\ I -Kj o 

since e, E md are positive, the matrix in {55) is Hurwitz for all Kp > 0, Kj > 0, and system l |55[ ) result 
asymptotically stable despite the unknown disturbance <po- 

The problem with implementing the regulator \55\ is that the resulting control signal is the average 
value of the actual control input rj, thus some procedure is required to synthesize a real- world control 
signal whose mean value satisfies the above conditions. In the SAF specific case this problem can be 
solved, consider the derivative of signal r\ a 

Jt^ = jMl T ^ dX (56) 

it can be rewritten on the left side as the difference over one period, while the right side is replaced with 
the derivative of r\ a expressed in {54); 

-fo (t) - 7] (t - T)] = -Kpi a + B = -Kpia - £K,Za (57) 

solving for rj (t) yields 

n{t) = -TKpt a {t)-sTK I z a {t) + n(t-T) (58) 

using {45), the derivative of the averaged square voltage error is actually measurable, thus the above for- 
mula is implementable. However, while the interconnection between voltage subsystem and controller 
is stable in sense of the averaged value, a further step is required. In the incremental implementation 
{58) there is no more an integral action, the control input history is kept in memory for one period, 
still the controller provides stability for the averaged voltage error z a - Consider now that for the phasor 



variables system, a stable steady-state guarantees that all the variables have a constant average value, 
while being allowed to oscillate freely. This property is desired for what concern the capacitor voltage 
and it is the main motivation for applying the averaging procedure, however implementation according 
to {58} can introduce undesired periodic oscillation in the control input 77, moreover oscillation will 
persist being remembered through the time delay term. In summary, while r\ a will approach the con- 
stant power loss value <po, the actual input T] might be any periodic signal with average value equal to 
<Po- Recalling that rj modifies the real power reference value x*,, any oscillation will result in a non-zero 
error for the power tracking controller. 

In order to avoid this situation the following term can be added to \5&\ 

d n {t) = T\{t-T)-T) a {t-T/2) (59) 

the reason of this modification is to cancel the oscillations stored in memory, by correcting the stored 
signal towards its own mean value r\ a (t — T/2). It is important to remark that the averaged value is not 
the actual mean value of its corresponding signal, the mean value s m of a signal s(t) is defined as 

1 rt+T/2 

Sm = - / s(n)dT (60) 

T Jt-T/2 

the above equation is identical to the zero-order phasor definition, except for a time shift ofT/2, For 
this reason the mean value of the stored signal tj (f — T) has been expressed as its time shifted average 
value, note that the mean value of this stored signal can be computed because also its "future" values 
are available. The additive term is a zero mean value signal, because it is obtained removing its 
DC-value to a periodic signal. Since the control input rj enters the averaged system ( |55) after being 
averaged itself, any modification having zero mean value will not affect the behavior of the averaged 
system dynamics. Hence the final implementation of control input together with \59\ is 

tj (t ) = -TK P ~z a - eTK lZ „ + r}a(t- 772) (61) 

this controller will not introduce undesired oscillation because it depends solely on averaged signals, 
whose simplified dynamics l |55[ > cannot give oscillations. 

4.2 Power tracking controller 

The voltage controller output reported in ( |61[ l can be replaced into the filter error power dynamics in 
\42\ , recalling also equation ( [54) , it turns out 

k = M{R,L)x- ~a+d(t)+f(e,za,e,ia,o) ( 62 > 

where 

d(t) = do+M(R,L)x*-x* +M(R,L)<po (63) 

is a T-periodic term composed by the sum of a constant term and sinusoids having known frequency, 
while 

f(z a ,e,Za,d,e) = TK p i a + eK I i a +K p i a {t-T/2)-S(t-T/2) 

(64) 

+M(R,L) [—TKpZ a - eKjZa -K p Za(t-T/2) + 0(t-T/2)]. 

The problem of forcing x in {62} clearly requires the ability of the control law to compensate for the 
signal d(t), perfect tracking cannot be achieved by a feedforward action since SAF parameters and 
d(t) are not fully known. To comply with uncertainties and provide robustness we propose an internal 



model-based controller. Each component of the vector d(t) can be seen as the output of the following 
linear system 

w i (t)=n Wi (t), Wi eTR 2N+1 

(63) 

dim(t) =r,w,-(f), i = d,q 

where T ; e r( 1x2/v+1 ) are suitably denned vectors and matrix O <E ]R(2tf+l)x(2tf+l) is defined as 
Q = blkdiag(Clj) with Oq = and 



o 



j(0„ 





, j=l,...,N 



(66) 



with the pairs (r,,n) observable. Defining <3> = bIkdiag(fi,Cl) and T = blkdiag(T d ,T q ), the following 
internal model-based controller can be designed 



« = r§ + Kx 



(67) 



where matrices Q and K need to be properly assigned. Once chosen u as in {67} and defined the internal 
model error variables as c, = t, — Lw, where w := [w^w^] 7 ^, the power subsystem closed-loop error 
dynamics can be rewritten as 



i=(M{R,L)--K)X--TZ+f(za,8,Za,6,e) 



(68) 



f = fl>f + 2x 



According to the general two time-scale averaging theory, the power tracking problem can be studied 
focusing on the boundary layer system, obtained by putting e = into the overall error dynamics, hence 

i becomes 



(69) 



by {47), d54b and z a = 0, 6 = 0, thus system 1 68 



% = (M(R,L) - \k)x- ~r| + f(za, 6,0,0,0) 



Now matrices K, Q need to be selected such that asymptotic stability is provided for the boundary 
layer system. Define two arbitrary Hurwitz matrices Fj, F q eR(2A' + l) x (2N+ 1), and two arbitrary 
vectors Gj, G q such that the pairs (Fj, Gj), (F q , G q ) are controllable, taking the controller matrices as 



K 



Q 



G d 
G 



K 



(70) 



with k[i, k q two arbitrary positive scalars, k a positive design parameter, and Ej, E q defined as non- 
singular solutions of the following Sylvester equations: 



Fd^d - E d Cl d — -GdT d 
FqEq F q Cl q — G q T q 



(71) 



asymptotic stability of the boundary layer system can be stated. In order to prove this claim let us define 
the vector 



R COmL 



(72) 



where I^i, I q \ denote the first element of vectors Tj, T q respectively and O2M is a zero raw vector 
having dimension 2N. Consider now the change of variables 



X = El-ER^ (0(t - T/2)-K p z a (t - Til)) + LGx 
where E = blkdiag(E c [,E q ), G = blkdiag{G,],Gq), in this coordinates system 1 69 1 results 

£= (M(R,L)- -K + TL~ l G)x- -TE^x 
% = F%- L(FG - GM(R,L) )x 



(73) 



(74) 



where F = blkdiag(F c i,F q ). Using standard linear system tools it can be verified that a value k exists, 
such that V k >k the state matrix of the system in the new coordinates is Hurwitz, hence asymptotic 
stability of the boundary layer system can be stated. 



4.3 Overall system stability 

Asymptotic stability has been stated for the boundary layer system and a linearized version of the 
averaged reduced voltage dynamics. Exploiting the main results of the two time-scale averaged systems 
theory, it can be proved that the two separately designed power and DC-bus voltage controllers, are able 
to provide practical stability for the complete system \A2\ , that is it's possible to claim that there exists 
a value £*, such that Ve < £*, k > k, I < I* the set 

{(jc,|) : x = 0, f = 0} x {(z, 6) : |f[ < /*, 8 = 0} 

is practically stable (Khalil 1996;) for the closed-loop trajectories of the complete error system. The 
fact that the proposed regulator is able to achieve the control objectives in a practical way means that 
the power vector x can be steered arbitrary close to the reference value while the averaged value r/ fl 
tends arbitrary close to the power loss term (fa. It's further to notice that the asymptotic tracking error 
can be arbitrary reduced by taking a smaller value for e, that is by increasing the capacitor value C. 



5. Simulation results 

Simulation tests have been performed in order to validate the proposed control solution. Two different 
scenarios have been adopted; first model {6]l has been implemented in MATLAB/Simulink and a load 
scenario with two harmonics at l(S>m and 13ft) m has been chosen. Then, in order to validate the con- 
troller performance in a situation closer to a real setup, the proposed continuous-time regulator has been 
discretized adopting a sampling frequency f s = 7 KHz, then the SAF converter components have been 
modeled by using Simulink/SimPowerSystems toolbox, and a suitable PWM technique with a carrier 
frequency equal to f s has been implemented. Finally a three phase diode bridge has been selected as 
nonlinear load scenario. 

The following system parameters has been set, according to the procedure illustrated in Section [3] 
C = 4400/iF, L = 33mH, R = 0.12O, while the DC-link voltage limits have been set to v,„ = 700V, 
vm = 900V. Ideal three-phase mains voltages with amplitude V,„ = 310V and frequency /,„ = 50Hz 
have been modeled. 

The internal-model based controller has been tuned to the load disturbances, according to the procedure 
described in |4.2| for what concern the simulations in time continuous domain. As regards the diode 
rectifier load scenario, the most relevant power disturbances, that is the 6 th and the 12 th load current 
harmonics expressed in the synchronous d — q reference frame (corresponding respectively to the 5'' 1 
and the 7''', and to the ll f * and the 13 th in the fixed reference frame), have been considered, then 
the IMC controller has been discretized according to the procedure reported in (Ro nchi et al.| [20031, 



thus the following matrices have been selected; Q = blkdiag(flo,ri(„rii2)> = T 9 = (1, 1,0, 1,0)^, 
K = diag(200,200) and Q = 10 3 diag{Q d ,Q q ), where Q d = Q q = (40.6,80.7,7.15,78.7, 17.6) r . For 
what concerns the voltage stabilizer described in |4.1| the following parameters have been selected 
K P = 0.3, Ki = 3.7. 

Consider now the performance obtained on the first simulation scenario, with ideal SAF model and the 
1 th , 13'' 1 disturbance harmonics; in Fig. |]the tracking error on both real and imaginary power variables 
is reported, as expected, asymptotic perfect tracking is achieved and the vector x is steered to the origin. 
This ideal behavior is confirmed by Fig. [7] [8j the two harmonics currents are totally canceled out by 
the filter currents, while a small current component oscillating at the first-order harmonic frequency and 
aligned to the corresponding voltage, arises on the line side due to the voltage controller action. In table 
|2]the harmonics compensation performance are summarized. 



Harmonic frequency [Hz] 


Una [A] 


Ha [A] 


Compensation percentage 


350 


0.0039 


10 


99.96% 


650 


0.0038 


10 


99.96% 



Table 2. Compensation performance for the two harmonics disturbance scenario. 

For what concerns the voltage controller, in order to validate the stability properties, a value quite far 
from the mean voltage reference value (vj n + vjL)/2 = 800 V has been chosen as initial condition for 
the capacitor voltage. As showed in Fig. [9] even though the average value is initialized at zero and needs 
one period before representing the actual voltage, the voltage controller reacts immediately, thanks to 
its dependance on the averaged derivative z a . Hence the voltage averaged error is successfully steered 
to zero, and the capacitor voltage is brought back to the middle of the safe interval, without exceeding 
the upper and lower bounds. The initial nonlinear behavior of the voltage error trajectories is originated 
by the neglected nonlinearities and also by the coupling term D(x), although it has been neglected 
due to two time-scale behavior hypothesis, it's excited by the internal model controller transient when 
harmonics compensation starts. 

As regards the second simulation scenario, carried out in discrete time domain and with a more detailed 
filter physical model, the power tracking performance are reported in Fig. [10] [TT] [T2] in this case 
the power error variables xj, x q are not exactly zero, due to the fact that AC/DC rectifier high order 
harmonics are not compensated by the internal model, furthermore the discretization effects have to 
be taken into account. However the load currents harmonics for which the controller has been tuned 
are strongly reduced at the line side as the currents magnitude spectrum reported in Fig. [12] shows. 
Analyzing the currents waveform in the time domain (Fig. 1 1 1[>, it can be verified that the mains currents 
are almost sinusoidal and perfectly aligned with the corresponding phase voltages, hence also the load 
imaginary power is almost totally compensated. The ripple introduced by the pulse with modulation can 
be noted on the filter current, it affects also the mains currents, however thanks to a correct inductance 
sizing, the high frequency distortion is properly bounded. Quantitative performance of the power- 
tracking controller obtained with this scenario are summarized in Tab. [5] 

The current component corresponding to the line frequency oscillation is slightly larger at the line side 
than at the load side, due to the additional active power drained to compensate for the filter losses. 
As regards the averaging voltage controller, a discrete time version has been implemented, while the 
same initial conditions of the first scenario have been reproduced. In Fig. [T3]the squared voltage error, 
it's averaged value and the actual capacitor voltage are reported, also in this case the objective relative 
to the voltage dynamics behavior is accomplished, similar considerations to those made for the previous 
scenario can be made. 



Harmonic frequency [Hz] 


ima [A] 


ha [A] 


Compensation percentage 


250 


0.03 


3.88 


99.2% 


350 


0.04 


1.91 


97.9% 


550 


0.03 


1.57 


98.1% 


650 


0.02 


1.08 


98.1% 



Table 3. Compensation performance for the diode bridge load scenario. 
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(a) Real power component tracking error, (b) Virtual power component tracking error. 
Fig. 6. Error variables xj, x q : two harmonics load scenario. 
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Fig. 7. Current and line voltage waveforms on phase a: two harmonics load scenario. 
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(a) Main current magnitude spectrum. 



(b) Load current magnitude spectrum. 



Fig. 8. FFT of the a-phase main current and of the corresponding load current: two harmonics load 
scenario. 



x 10 




0.2 0.4 

time [s] 



0.2 0.4 

time [s] 



0.6 



(a) Square capacitor voltage error and com- (b) Actual capacitor voltage value, 
puted average value (bold). 



Fig. 9. Voltage controller performance: two harmonics load scenario. 
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(a) Real power component tracking error, (b) Virtual power component tracking error. 
Fig. 10. Error variables xj, x q : diode bridge load scenario. 
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(a) Main current magnitude spectrum. (b) Load current magnitude spectrum. 

Fig. 12. FFT of the a-phase main current and of the corresponding load current: diode bridge load 
scenario. 
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Fig. 13. Voltage controller performance: diode bridge load scenario. 



6. Conclusions 

In this chapter a nonlinear robust control solution for a shunt active filter has been proposed, the focus 
has been firstly put on the hardware components design issue, providing a suitable algorithm, based on 
the structural system properties, which gives guarantees on the feasibility of the control problem and 
allows to obtain a crucial time-scale separation between the power and voltage dynamics. Then ex- 
ploiting nonlinear systems analysis well established tools, such as averaging and singular perturbation 
theory, an averaging capacitor voltage controller and a power tracking controller based on the internal 
model principle, have been presented. The former exploits the insight that, regulating the averaged 
voltage value, makes it possible to ignore the necessary oscillations for a proper filter operation, and 
improves the voltage dynamics behavior. The second is chosen in order to ensure asymptotic tracking 
of undesired load current components, providing also robustness with respect to disturbances and model 
uncertainties. 

Saturation issues have not been explicitly addressed in this work, owing to space limitation, however it 
is of utmost importance to deal with these phenomena for an actual industrial implementation with sta- 
bility and performance guarantees. Some solutions, for the SAF specific case, have been proposed (see 
|Cavini et al.| ( |2004") ,|Cavini et al. (2004)), however this is still an open research topic. Future effort will 
thus be devoted to improve the filter performance under control input saturation, analyzing the problem 
in the context of modern anti-windup approaches, hence providing a rigorous characterization of the 
system under saturation constraints. Moreover discretization issues relative to the nonlinear controller 
here discussed will be further analyzed, in order to improve the discrete-time controller performance 
with respect to that obtained applying standard discretization techniques. 
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